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ABSTRACT 

We explore trends in galaxy properties with Mpc-scale structures using catalogues of 
environment and large scale structure from the Galaxy And Mass Assembly (GAMA) 
survey. Existing GAMA catalogues of large scale structure, group and pair membership 
allow us to construct galaxy stellar mass functions for different environmental types. To 
avoid simply extracting the known underlying correlations between galaxy properties 
and stellar mass, we create a mass matched sample of galaxies with stellar masses 
between 9.5 ^ logM»//i“^M 0 ^ 11 for each environmental population. Using these 
samples, we show that mass normalised galaxies in different large scale environments 
have similar energy outputs, u — r colours, luminosities, and morphologies. Extending 
our analysis to group and pair environments, we show galaxies that are not in groups 
or pairs exhibit similar characteristics to each other regardless of broader environment. 
For our mass controlled sample, we fail to see a strong dependence of Sersic index or 
galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. 
Repeating our analysis for galaxies that have not been mass controlled introduces and 
amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, 
indicating that stellar mass is the most important predictor of the galaxy properties 
we examine, as opposed to environmental classifications. 

Key words: galaxies: evolution - galaxies: clusters: general - galaxies: luminosity 
function, mass function - galaxies: stellar content - large-scale structure of Universe 


1 INTRODUCTION 


No galaxy is an island. The close surroundings of a galaxy 
(i.e. within a few tens of Mpc) play an important role in 
shaping its evolution (e.g. Avila-Reese et aT]|2005 Blanton 


et al.|2005 Hahn et al.|2067 Fakhouri k. Ma|2009 l. Stellar 

populations are particularly sensitive to environment: the 
proximity of nearby or merging galaxies has been shown 


to trigger gas collapse within galaxies, which leads to an 


increased rate of star formation (Toomre & Toomre 


1972 


Quinn et al. 

|1993| |Lewis et al.||2002| |Gomez et al. 2003 

Porter et al. 

20081. Other observational signatures of envi- 

ronmental effects include differences in colour (|Kreckel et al. 

20121, stellar mass (Chabrier 20031, morphology (Dressier 


1980 Butcher fc Oemler||1984 Dressier et al. 19971, the In- 
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2 M. Alpaslan et al. 


minosity function (Croton et al. 2005 McNaught-Roberts 


et al. 20141, and disrupting the gas content of a galaxy 
( Beygu et al.|20l3 Benitez-Llambay et al.|2013 |. 

Beyond the close surroundings of a galaxy, there is also 
ample evidence to show that, for galaxies in a cluster, the 
cluster environment has some influence on its properties. 
Galaxies do appear to be more red in denser environments 


(e.g. Baldry et al. 2006 Bamford et al. 2009 Peng et al. 


2010 ). On the other hand, star formation appears to be sup¬ 
pressed in galaxies that are infalling into clusters (|Lewis| 


et al. 2002 


Balogh 20071; though more recent work by 


Brough et al. 


(20131 fails to detect significant trends in star 


formation rates with environment. Similarly, recent work by 


Wijesinghe et al. (20121 determines that star formation rates 


of morphologically classified galaxies depend, to first order, 
on stellar mass and not environment. Another well known 
phenomenon is the tendency for galaxies in different group 


environments to have distinct luminosity functions (e.g. Nor- 
berg^^t^lJ[2002jJDePropris^^et^ay[2003|J^otra^etalJ[2005|^ 


Zehavi et al. 2005j^ Robotham et al. |2006[ 2010^ JZehwr 


et aQ 201^ Maaaki_e^^^ 2013| IMcNaught-Roberts et al. 

llEardlV 


2014[|Eardley et al.|2014p' 

Galaxies, however, do not just live amongst their im¬ 
mediate neighbours; instead, they are embedded within the 
complex large scale structure of the Universe. This struc¬ 
ture is characterised by long, linear filaments of galaxies that 
stretch between highly dense clusters containing a vast num¬ 
ber of galaxies. In between filaments and clusters lie voids, 
which are largely empty regions of space that contain a mere 
handful of galaxies per unit volume. Does this so-called ‘Gos- 
mic Web’ have any noticeable impact on galaxies? In other 
words, is a galaxy in a filament distinctly and systematically 
different to that inside a void? 

Research in understanding the Gosmic Web is actively 
ongoing. A number of algorithms and methods exist to de¬ 
tect and classify large scale structure (e.g. 


Doroshkevich 


et al. 

20041 Pimbblet 20051 Golbergl 

20071 IForero-Romero 

et al.| 

20091 IMurphy et al. 1120111 ISmit 

1 et al. 1120121 IGautun 


et al.|2012 [Alpaslan et al.|2014a [Eardley et al.|2014| ). The^ 

algorithms primarily focus on using positional (and where 
possible, velocity) information in the distribution of galax¬ 
ies to trace and identify not only filamentary structures, but 
also voids and clusters (often referred to as ‘knots’). These 
sophisticated algorithms, together with current generation 
redshift surveys make it possible to compare and correlate 
the properties of galaxies found in a range of environments. 
The Galaxy and Mass Assembly (GAMAQ survey 


(Driver et al. 2009 2011| Liske et al. 2014) is a joint 


European-Australian project which combines multiwave¬ 
length photometric data from a large number of ground and 
space-based programs with a highly complete spectroscopic 
campaign conducted using the AAOmega spectrograph on 
the Anglo-Australian Telescope at Siding Spring Observa¬ 
tory, Australia. At the time of writing, GAMA has access 
to data in the following bands (see Driver et al. in prep): 
Galaxy Evolution Explorer (GALEX; FUV, NUV), Sloan 
Digital Sky Survey (SDSS; ugriz), UKIRT Infrared Deep 
Sky Survey (UKIDSS; YJHK), VISTA Kilo-degree Infrared 
Galaxy survey (VIKING; ZYJHK), Wide-Field Infrared 


http: / / www.gama-survey.org/ 


Survey Explorer (WISE; Wl-4; MIR), and Herschel ATLAS 
(H-ATLAS; PACS, SPIRE; FIR) and Giant Metrewave Ra¬ 
dio Telescope and shortly Australian SKA Pathfinder (ra¬ 
dio). The spectroscopic component of the survey contains 
~250,000 spectra for galaxies out to r = 19.8 mag across 5 
fields: centred at q = 9h, S = 0.5 deg (G09), a = 12h, S = 
-0.5 deg (G12) and a = 14h30m, S = 0.5 deg (G15), a = 
2h, 5 = -8.125 deg (G02) and a = 23h and <5 = -32.5 deg 
(G23). The three equatorial fields (GOO, G12 and G15) are 
12 X 5 degrees each, and the two southern fields (G02 and 
G23) are respectively 8.6 x 2.5 and 12 x 5 degrees. 

One of the principal advantages of GAMA for studies of 
intergalactic structure is the high spectroscopic target den¬ 
sity; which averages to 1,050 galaxies per square degree, and 
a spectroscopic completeness of 98.42% in the three equato¬ 
rial fields, with associated velocity uncertainties of « 50 
kms“^ (Liske et al. 2014). This is achieved via an observing 
strategy which revisits the same patch of sky, on average, 
10.3 to 10.9 times (Liske et al. in prep), with different fi¬ 
bre configurations on the 2dF instrument, in order to ob¬ 
serve every possible target ( Robotham et ar]|2010a |. These 
highly complete and dense data have led to the creation 


of th e GAMA Galaxy Group Gatalogue (Robotham et al. 
G^G) and the GAMA Large Scale Structure Gata- 


2011 


logue ([Alpaslan et al. 2014a| GLSSG), both of which will 


be discussed in greater detail in Section These data sets 
allow us to tackle outstanding questions in galaxy evolu¬ 
tion from an empirical perspective, and will help in placing 
robust constraints on how the large scale structure of the 
Universe affects galaxy evolution; particularly morphology, 
colour, and stellar populations. Ultimately a detailed study 
of the gas content will also be required, possible via projects 
such as the ASKAP DINGO survey ( Meyerf2009 1. 


In this work, we seek to make some initial advances into 
understanding to what extent (if at all) large scale structure 
has an impact on the evolution of galaxies. Given the broad¬ 
ness of this topic, we choose to focus on a few, relatively well 
understood properties of galaxies: colour, brightness, mor¬ 
phological properties (Sersic index, ellipticity, effective ra¬ 
dius, and visually identified structure), gas metallicity, and 
total spectral energy density. For each of these properties, 
we seek to find any possible trends related to large scale 
structure. In other words, do galaxies that exist in different 
types of large scale structure inhabit different parts of the 
parameter space of galaxy properties? How do these trends 
compare to those of galaxies in groups, or pairs? 

In Section we give further detail on the G^G and the 
GLSSG, as well as other GAMA catalogues of the various 
galaxy properties that were used in this investigation. In Sec¬ 
tion 3, we display galaxy stellar mass functions for GAMA 
galaxies in different environments. Section[^begins by defin¬ 
ing a mass normalised sample of galaxies, which we then 
use to look for trends in galaxy properties within different 
types of environment. Finally, in Section 5 we briefly ex¬ 
amine consequences of our mass normalisation process and 
repeat our analysis for different mass ranges. In Section 6 
we discuss these results and summarise them. Throughout 


this paper, consistent with the cosmology used in 

RobothamI 

et al. (20111; Merson et al. 

(20131; Alpaslan et al. ( 

2014a 

), we 


All of the analysis in this work has been done using the R 
programming language (R Gore Team|[2014K 
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Galaxy properties in large scale environments 3 


2 DATA 


The large scale catalogue we use is the GLSSC (Alpaslan 
et al.|[^14a| FilamentFindingvOl), an observationally mo¬ 
tivated catalogue of large scale structure. It classifies galax¬ 
ies in the three equatorial GAMA fields as belonging to 
hlaments or voids. The catalogue also introduces a third, 
interstitial structure of galaxies dubbed ‘tendrils,’ which 
typically consist of delicate (but coherent) strings of 6-10 
galaxies that emerge from filaments and terminate either in 
voids, or in other filaments. Tendrils are discussed in greater 


detail in Alpaslan et al. (2014b I. The catalogue contains 


45,542 galaxies with Mr ^ —19.77 + 51og/i mag across all 
three equatorial fields, each assigned to a filament, tendril, 
or void. The full sample of galaxies is volume limited to 
Mr SS —19.77 + 51og/i mag and extends out to z = 0.213. 
Gonsidering the GAMA selection cut of r = 19.8 mag, this 
redshift and its associated absolute magnitude cut are cho¬ 
sen such that they maximise the number of galaxies in the 
GLSSG. The catalogue also identifies 643 filaments, com¬ 
posed of an average of 8 groups and spanning u p to 100 
h”'^ Mpc. The appendix of Alpaslan et al. (2014aI goes into 


further detail on the actual composition of the catalogue. 

At its heart, the algorithm that generates the GLSSC 
relies on a slightly modified minimal spanning tree (MST) 


Golberg (2007 

and uses a multiple pass approach, similar to 

Murphy et al. 

(20111. Initially, a volume limited sample of 


ning tree to identify filaments; galaxies within a distance r of 
each filament are then associated with that filament. These 
galaxies, within and outside of groups, are together referred 
to as filament galaxies. The remaining galaxies (not in fil¬ 
aments) are processed through a second minimal spanning 
tree to identify tendrils and tendril galaxies; and any galaxy 
beyond a distance q from a tendril is considered to be a 
very isolated galaxy, which we call a void galaxy. Note that 
this approach makes no attempt to physically identify voids. 
The values for r and q are chosen such that the two point 
correlation function ^ 2 (r) of void galaxies is minimised i.e. 
our fiducial assumption is that voids should have a minimal 
amount of coherent structure due to the low net accelera¬ 
tions experienced over the lifetime of the Universe. A third 
parameter, b represents the maximum allowed distance be¬ 
tween groups in order for them to be considered to be in the 
same filament; this is chosen such that at least 90% of groups 
with L* ^ 10 ^^Lq are located in hlaments. The GLSSC is 
generated with b = 5.75h~^ Mpc; r = 4.12h~^ Mpc; and 
q — 4.56/i“^ Mpc. 


2.1 Additional GAMA catalogues 

Here we briehy introduce each of the other catalogues that 
are used for the analysis in this work. Where possible, we ref¬ 
erence the publication that accompanies the catalogue, and 
suggest that the more curious reader refers to these for addi¬ 
tional detail. Where no reference is given, the catalogue has 
been provided by one of the authors of this paper. We con¬ 
sider only data for galaxies in the three equatorial GAMA 
helds. 


2.1.1 Galaxy spectra 


This catalogue (Liske et al. in prep, SpecCatv25), contains 
spectra for all galaxies observed using the 2dF instrument 
at the Anglo-Australian Telescope in Australia, for GAMA. 
Spectra obtained at the AAT were redshifted by observers 
at the telescope, as well as via an automated algorithm Au- 
TOZ (see Baldry et al.|2014 for details on Autoz, and Liske 
et al. in prep for details on the spectroscopic campaign). On 
average, redshift measurements for the spectra have an asso¬ 
ciated error a of approximately 27 km s“^; and only 0.2% of 
redshifts are considered to be incorrect (Baldry et al.|[2014 


Liske et al. in prep). Some spectra from previous surveys 
(e.g. SDSS) are also included in this catalogue. These spec¬ 
tra (and the redshifts derived from them) form the most 
important data asset for the catalogues of environment dis¬ 
cussed below. 


2.1.2 The GAMA Group Catalogue (G^C) 


The G'^’C (Robotham et al. 


2011 


GroupFindingv07) is 


group catalogue that has been put together using a slightly 
modified friends-of-friends algorithm that considers galaxies 
to be in a group if they are grouped both along the line of 
sight, as well as when projected onto the sky. This success¬ 
fully accounts for redshift space distorsions caused by the 
peculiar velocities of galaxies in groups. The catalogue con¬ 
tains 23,838 groups with 2 or more members out to r = 19.8 
mag; these groups contain approximately 40% of all galax¬ 
ies in GAMA. This catalogue also provides galaxy pairing 
information, which is also used in this work and studied in 
[Robotham et al.] ( |2014[ ). Parameters for the groupfinding al¬ 
gorithm is calibrated by being run on a series of GAMA 


mock galaxy catalogues, described in Merson et al. (20131, 


and the optimal linking lengths have been verified as theoret¬ 


ically optimal in recent work by Duarte & Mamon (20141. 


The mock catalogues are designed to mimic the geometry 
of the fields observed by GAMA, as well as replicate the 
galaxy luminosity function of the survey. They are assem¬ 
bled by populating haloes taken from the Millennium Sim¬ 
ulation ( Springel et al.|[2005 l using the Bower et al. (20061 
Galform semi-analytic galaxy formation model. 


2.1.3 Aperture matched photometry 


Aperture matched photometry for GAMA (Liske et al. in 
prep, ApMatchedCatvOS) is measured by running Source Ex¬ 


tractor (Bertin & Arnouts 1996 Berlin,E. 20111 on imag¬ 


ing data from the SDSS and UKIDSS (for ugriz bands and 
Y JHK bands respectively) in matched aperture mode. This 
allows us to place apertures on sources in the r-band (as this 
has the highest imaging quality and defines the main sur¬ 
vey spectroscopic sample selection) and measure fluxes in 
all bands. The imaging data have been preprocessed into a 
series of 27 (9 bands x 3 fields) large mosaics using SWarp 
( Bertin et al. [20(12 ( and have been normalised to a common 
zeropoint and convolved to a common PSF of 2 arcseconds 
(see Hill et al.|2011 and Driver at al. in prep). 
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2 . 1.4 Stellar masses 


The GAMA stellar masses (Taylor et al. 2011 


StellarMassesvl6) are estimated using photometry in 
the restframe ugriz bands from the aperture matched 
catalogue ApMatchedCatv05, to which a series of synthetic 
spectra are fit. The spectra are designed to incorporate 
an exponentially decaying star formation history, using 


Bruzual & Chariot (20031 models with a Chabrier (20031 


IMF and the Calzetti et al. (20001 dust obscuration law. 


The stellar masses are determined by integrals that are 
weighted to the probability of each fit. This catalogue 
also provides dust and extinction corrected apparent and 
absolute magnitudes, as well as restframe and extinction 
(Galactic and self-attenuation) corrected u — r and g — i 
colours, which are also used in this work. The extinction 
correctios are estimated during the sythentic spectral hts 
and have associated uncertainties of « 0.15 mag (Taylor 


et al. 20141, but this is likely to be an overestimation as 
the random errors on the rest-frame colours are on the 


order of 0.05 mag. As explained in Taylor et al. (20141, the 


amount of dust from which A„ is derived for each galaxy 
cannot be negative, leading to a possible overstimation of 
Ay for dust-poor galaxies. This will translate into an error 
in [g — i) colours of up to 0.1 mag; however, these errors 
have no discernible impact on our results. 


2.1.5 Sersic photometry and morphology 


This catalogue provides single component Sersic hts 
of galaxies, computed (SersicPhotometryvO?) using the 
SIGMA software which wraps around various astronomy 
packages including GALFIT3 ( [Kelvin et al.|201^|Peng et al.| 
2011). Along with providing a Sersic index for each 


galaxy, this catalogue also provides some other morpholog¬ 
ical parameters such as effective radius Re and ellipticity, 
dehned as 1 — b/a, where a is the ratio of the semi-major 
and b the minor axes of the galaxy. An ellipticity value of 
0 being given to galaxies whose semi-major and semi-minor 
axes are equal (i.e. round galaxies). 


2.1.6 Emission line measurements and metallicities 


This catalogue measures fluxes and equivalent widths for 
spectra for all GAMA galaxies. Measurements are made 
for 11 different emission features, including Ha, H/3, Olll 
A5007, and Nil A6583, using single and double Gaussian hts 
to continuum-substracted spectra. This is done in a manner 
similar to that described in [Gunawardhana et al.j ( |2013| ). 
Gas metallicities of galaxies are then calculated using the 
03N2 index. This is dehned as: 


03N2 = logj^o 


/ [Oiii] A5007/H/3\ 
[Nil] A6583/Ha ) 


( 1 ) 


from which it is possible to calculate the O/H metallicity 
indicator as [12 -|- log(0/H)] = 8.73 — 0.32 x 03N2 as per 
the prescription of jPettini fc Pagel] (2004). Line strengths 
are emission measurements and have been corrected for dust 
attenuation, and have been calibrated to match metallicities 


given by Tremonti et al. (20041, as per the prescription given 
inf 


Lara-Lopez et al. (2013). 


2.1.7 GALEX photometry 

This catalogue provides GALEX NUV and FUV photome¬ 
try for the GAMA II equatorial survey regions. The GALEX 
ultraviolet catalogue (GalexPhotometryv02; Andrae et al. in 
prep) is a combination of archival data and pointed obser¬ 
vations on equatorial GAMA helds. The archival data have 
been used to extend the ultraviolet coverage of the GAMA 
regions as much as possible. 


2.1.8 WISE photometry 


Photometry for the four mid-infrared WISE bands is pro¬ 
vided in WisePhotometryvOl for GAMA galaxies in all three 
equatorial helds that are matched to WISE observations. 
The catalogue does not include sources that have no GAMA 
match within 3 arcseconds of a WISE source (see [Gluverj 
et al. 2014 for details). The photometry for each source is 
from the All-Sky WISE Data Release with standard aper¬ 
ture measurements for unresolved sources and isophotal pho¬ 
tometry for resolved sources. Prohle-ht measurements are 
also provided. 


2.1.9 Multi-band photometry 


The process of amalgamating data from multiple sources 
into a single photometry catalogue (20BandPhotomv02) is de¬ 
scribed in Driver et al. (in prep) and builds on the work de¬ 
scribed in Hill et al. ( 2011| ); here we provide a brief summary. 
The GALEX (FUV, NUV), aperture matched photometry 
{u-K) and WISE {Wl-4) catalogues are combined to PAGS 
and SPIRE photometry derived from the PAGS/SPIRE 
maps by measuring the hux in the appropriate optically de¬ 
hned aperture convolved with the appropriate PAGS/SPIRE 
psfQ This combination is done by exact name matching 
to unique identihers that are given to each matching galaxy 
using TOPCAT ( Taylor|2005 1. The combined catalogue is cor¬ 
rected for Galactic extintion using E{B — V) values provided 
by Schlegel et al. (19981 and the coefficients listed in Liske et 


al. (in prep). The hnal product is a catalogue that contains 
invidual hux measurements in Janksys across all bands for 
each individual GAMA galaxy, with dummy hux values in¬ 
cluded if that galaxy has not been surveyed in a particular 
band. This data is used only in Section [4.3| 


2.1.10 Visual morphology 

For a sub-sample of galaxies with 0.002 2 ^ 0.1, r < 19.8 

mag, this catalogue (VisualMorphologyv02) provides visual 
classihcations of galaxy morphology as being either elliptical 
or not-elliptical. The classihcation was performed on postage 
stamps generated from three colour giH band images with 
arctan scaling from the SDSS and VIKING or UKIDSS LAS 
data. The postage stamp sizes were set to a constant value 
of 30 kpc X 30 kpc, except if the implied size is greater than 
100 pixels, in which case it was set to that value. 

^ Where objects overlap care is taken to divide the flux between 
the two objects following the prescription outlined in Appendix A 
of [Bourne et al.] l |2012| | (using the aperture matched photometry 
catalogue as the input catalogue to deblend data), see also Driver 
et al. in prep. 
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2.1.11 Surface density 


The nearest neighbour surface density, E5 (Brough et al. 
2013[ EnvironmentMeasuresv5), is calculated for all GAMA 
galaxies with reliable redshifts (nQ > 2; Driver et al.|20lT |. 
The surface density is dehned using the projected comoving 
distance to the fifth nearest neighbour (d5) with ilOOOkms^ 
within a volume-limited density-defining population: E5 = 
5/7rd|. The density-dehning population has absolute SDSS 
petrosian magnitudes < Mr,iimitQz, k-corrected to z = 


0 following Loveday et al. (20111, where Mryimit = 20.0 mag 
and Q defines the expected evolution of Mr as a function of 
redshift (Q = 0.78; Loveday et al. in prep). Densities are 
then corrected for the survey r-band redshift completeness 
as E5 = E5,rau) X 1/completeness. Galaxies where the nearest 
survey edge is closer than the fifth nearest neighbour are 
flagged as upper limits. 


2.2 Sample selection 


Galaxy properties are well known to correlate with stellar 
mass. To ensure we recover variations due to environment 
rather than these known stellar mass dependencies, it is vital 
that we hrst remove the mass dependency from our sample. 
We do this by forcing our samples to be mass complete, 
using stellar mass data from the StellarMassesvl6 cata¬ 
logue. We choose the limits a 0.1 and M* ^ 10®'® 
such that we can use the maximal number of galaxies from 
all catalogues considered in this work, each with its own 
redshift limits. Our sample is shown in Figure where it 
is highlighted according to each galaxy’s restframe g — i 
colour, and the red lines mark our limits. In |Alpaslan et ah] 
1 2014a I it is established that the GLSSC is mass complete 
only to logM*/h“^M0 ^ 10.61 (but does include galaxies 
of lower mass than this). To account for the discrepancy 
between this mass limit and our chosen sample, we take 
all galaxies with \ogMt./h~^M q ^ 9.5 with z 0.1, and 
Mr ^ —19.77 -I- 51og/i mag and classify them as belonging 
to the same type of large scale structure (i.e. filament, ten¬ 
dril, or void) as their nearest neighbour using the GLSSC. 
Our low mass large scale structure sample contains 7195 
volume limited galaxies. 

In [Alpaslan et al.| | |2014b| we show that void galax¬ 
ies in the GLSSC are isolated and not associated with any 
structures by making use of the line correlation function. 
Introduced in [Obreschkow et al.| ( [2M^ , the line correlation 
function l{r) measures the amount of linearity in structures 
of length r in a way that is analogous to how the two-point 
correlation function measures clustering on length scales of 
r. The line correlation function takes the phase factor in¬ 
formation of the galaxy density held (i.e. e(fc) = 5{k)/\5{k)\, 
where 5(k) is the Fourier transform (FT) of a galaxy density 
held 5(f)) and calculates the normalised three-point corre¬ 
lation of the inverse FT of i{k) for three points on a straight 
line separated by a distance r. For an illustration of how 
the line correlation function can distinguish between ran¬ 
dom, spherical, and linear overdensities, we refer the reader 


to Figure 3 in Obreschkow et al. (20131. 


We repeat this calculation for the galaxies in our low- 
mass GLSSC sample in order to conhrm that the process 
of nearest galaxy association described above has not re¬ 
sulted in the association of galaxies in structures with void 


CO 


0 0.1 0.2 0.3 0.4 0.5 0.6 

Z 



Figure 1. Stellar mass as a function of redshift for all GAMA 
galaxies. Galaxies with z ^ 0.1 and AL* (S lO® ®/i“^AL 0 are 
coloured by their restframe g — i colour. Redshift and mass limits 
are shown as red lines. These limits are chosen as they are shared 
by all the catalogues we consider in this work. 


galaxies. As was done in Alpaslan et al. (2014bI, we mea¬ 
sure the filamentarity of each galaxy population by calcu¬ 
lating its so-called excess line correlation with respect to 
a random points set. The excess line correlation is defined 
as Al(r) = [Z(r) — loir)]^/f, where l{r) is the line correla¬ 
tion of galaxies in the GLSSC placed into a box of side 150 
h“^ Mpc, Zo(r) is the line correlation of an equal number of 
galaxies placed randomly into a cone identical to the GAMA 
geometry within this box, and / is the volume fraction of 
the GAMA fields within this box. Including / ensures that 
Al{r) is approximately independent of the volume of the box 
(see Section 3.4 in [Obreschkow et al.|2013 |. 

In Figure we show the excess line correlation of galax¬ 
ies in the low mass GLSSC sample averaged across all three 
equatorial GAMA fields. Notably, filament galaxies continue 
to have the highest Al(r) signal, while void galaxies show no 
filamentarity. Errors shown are the variance of Al(r) across 
all three GAMA fields, with the scatter due to the smaller 
number of points in this sample. The lack of filamentarity 
in the void galaxy sample confirms that we have not inad¬ 
vertently associated galaxies within structures with the void 
galaxy population. Figure conhrms the linearity (or lack 
thereof) present in the three large scale structure classifica¬ 
tions. 

We further test our method by taking the GLSSC and 
randomly removing 50% of the galaxies in each environment 
(hlament, tendril, and void). We then apply the reclassi- 
hcation scheme described in the beginning of this section 
to these removed galaxies, and compare that classification 
to what the GLSSC classihcation is. This random rejection 
process is realised 100 times, and the average recovery rate 
with the correct classification is 90% for all three categories 
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Figure 2. The excess line correlation Al(r) of galaxies in fila¬ 
ments, tendrils, voids, and all galaxies shown in blue, green, red, 
and black respectively. We compute Al{r) individually for each 
population of galaxies in all three GAMA fields, and show the 
averaged results in this figure. The errors shown are the standard 
deviation in the excess line correlation across the three GAMA 
fields. Most notably, we still retain the 0 signal for void galaxies 
across large distances. Readers making a direct comparison be¬ 
tween this figure and Figure 2 in [Alpaslan et al.| | |2014b| | must 
note that the line correlation calculation method has since been 
updated (see Wolstenhulme et ah, submitted for details). This 
is mainly reflected by a change in the vertical normalization of 
Al(r). 


of large scale structure (with a standard deviation sub-1%, 
indicating that this process is robust to randomness). This 
high recovery rate is a strong indicator that the reclassifica¬ 
tion scheme does not create large biases in our large scale 
structure sample. The « 10% of galaxies that are misclassi- 
fied are usually those whose neighbours were also randomly 
removed, resulting in the nearest structure being of another 
type. Their spatial distribution is randomised within each 
realisation, and the mean distance between a reclassified 
galaxy and its nearest neighbour is 1.97 h~^ Mpc for this 
test case, and approximately 10% greater for the reclassifi¬ 
cation scheme used to generate the low mass sample. 


3 GALAXY STELLAR MASS FUNCTIONS BY 
ENVIRONMENT 

In order to broaden our investigations beyond the simple fil¬ 
ament, tendril and void sub-samples we also look to now 
incorporate some of the additional information available 
within the GAMA database, i.e., pairings and groups. Fig- 
ure[^ shows t he galaxy stellar mass functions (GSMFs; see 
Baldry et al. 2012 Kelvin et al. 20141 of galaxies in a va¬ 


riety of different environments, including those in the low 
mass GLSSG sample. The environmental types for which we 


calculate GSMFs in Figure are; filaments, tendrils, voids 
(as defined in the GLSSC), galaxies in groups with high 
mass halos, mid mass halos and low mass halos (defined as 
Mh > 10^"^M0, 10^® and Mh ^ 10^®Mg 

respectively), galaxies not in groups, galaxies in pairs, and 
galaxies not in pairs (any two galaxies with a physical pro¬ 
jected separation of 100 h~^ kpc and velocity separation of 
1000 kms“^ are considered to be in a pair). Following the 


prescription of Kelvin et al. (20141, we calculate the GSMF 


in the form of a Schechter function: 

<E>(M) AM = ln(10)(;A*10'^“"*^^“+^’ exp (^-10^"""*^^ dM(2) 

where M = logj^Q M. The fit parameters for the stellar 
mass functions are given in Table The errors in the 
parameters are estimated via jackknife resampling, with 

where x is the best fitting pa¬ 
rameter and Xj is the best fitting parameter as given from a 
jackknife resampled version of the data. The jackknifing is 
done in such a way as to split the data into A = 10 equal 
sized bins. Across all populations, the slope at the low mass 
end a is statistically consistent with being constant. The rise 


at the low mass end of the overall GAMA GSMF (Baldry 
et al.|20T2 1 must therefore be caused by new low-mass pop¬ 


ulations being sampled, and not due to any fluctuations in 
the numbers of other populations below M*. 


ARE GALAXY PROPERTIES INFLUENCED 
BY LARGE SCALE STRUCTURE? 


4.1 Mass normalisation 

From Figure it is evident that galaxies in different envi¬ 
ronments have vastly different stellar mass functions. Given 
that stellar mass most directly influences galaxy evolution 
(e.g. Oesch et al.||2010 Robotham et al.]|2013 l, performing 
our analysis on our galaxy sample as is is more likely to draw 
out trends in stellar mass; or conversely, blur out any trends 
caused by environment. 

This mass segregation can be accounted for by resam¬ 
pling each environmental population of galaxies such that 
their galaxy stellar mass functions all match. In other words, 
we resample ^ch population such that within the interval 
[M, , M* -I- dM*] is equal across all galaxy populations. We 
split each population into bins of size 0.05 dex of stellar mass 
between 9.5 ^ logM,/h~^ Mq ^ 11, and in each bin, ran¬ 
domly resample the galaxies in each population (except for 
the voids) so that their GSMF in that bin matches that of 
the void population in that same bin. Galaxies whose mass is 
greater than 10^^ log Mt/h~^ Mq are discarded. The results 
of this process can be seen in Figure where we plot the 
PDFs of stellar masses of galaxies in different environments. 
We choose to match all populations to the GSMF of the void 
galaxy population, as this is the most extreme (and because 
we cannot ‘upscale’ that GSMF to the others). The stellar 
mass normalised samples are therefore random subsets of 
each galaxy sample that have the same mass distribution as 
the void sample. From here on, whenever we refer to a partic¬ 
ular population of galaxies, we always refer to a sub-sample 
with mass matching; meaning that we are discussing galaxies 
whose stellar mass is between 9.5 ^ logM*/fe“^Mg ^ 11. 
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Table 1. Parameters for the single Schechter function GSMF fits shown in Figure]^ The columns are the knee in the Schechter function 
(M*), the slope at the low mass end (a), the normalisation {<p*; in units of dex“^ Mpc“®), and the goodness of fit jv. The fractional 
integrated stellar mass from the fits is also shown, for subdivisions of environment (marked by the horizontal lines in the table). Errors 
are estimated from jackknife resampling. 



logM*/M 0 

OL 

0 */lO"® 


f ^(M) dM 

All 

10.82 ± 0.02 

-0.97 ± 0.02 

2.00 ±0.13 

1.16 

100 % 

Filaments 

10.88 ± 0.02 

-0.99 ±0.02 

1.52 ± 0.07 

1.2 

57% 

Tendrils 

10.66 ±0.03 

-0.86 ±0.05 

1.55 ± 0.01 

1.09 

39% 

Voids 

10.31 ± 0.55 

-1.24 ± 0.62 

0.15 ± 0.19 

2.31 

4% 

HlMass 

11.16 ± 0.10 

-1.06 ±0.07 

0.27 ± 0.06 

1.25 

10 % 

MidMass 

10.93 ±0.025 

-0.87 ± 0.02 

1.12 ± 0.05 

1.59 

24% 

LoMass 

10.86 ± 0.02 

-0.87 ± 0.02 

1.74 ±0.06 

1.16 

36% 

Ungrpd 

10.64 ± 0.04 

-1.07 ± 0.03 

1.35 ± 0.11 

1.60 

30% 

inPair 

10.89 ±0.02 

-0.87 ± 0.03 

1.22 ± 0.06 

1.28 

30% 

notinPair 

10.80 ± 0.02 

- 1.02 ± 0.02 

2.50 ±0.12 

1.01 

70% 


Figure shows the relative fraction of galaxies in their envi¬ 
ronments, split by large scale structure, group membership, 
and pair membership. 

For the remainder of this paper, we will combine our 
mass normalised environmental populations of galaxies with 
the various catalogues introduced in the preceding section. 
The sample we use is the galaxies that are common to all 
catalogues and are within our specified mass and redshift 
ranges. When mass normalised, this sample contains 28,781 
galaxies. The sample size increases to 35,850 galaxies when 
the mass normalisation is removed. Note that all environ¬ 
mental categories of galaxies are mass matched indepen¬ 
dently. 


4.2 Colour, brightness, and morphology 


Recent work studying the properties of void galaxies has 
shown them to be predominantly bluer and fainter with re¬ 
spect to their counterparts in higher density regions (e.g. 
Kreckel et al. 20121. Conversely, work by Eardley et al. 


(20141 which uses a tidal tensor description to classify 


GAMA galaxies into knots, filaments, sheets, and voids finds 
that the luminosity function of galaxies is independent from 
their large scale structure classification. In this work, we 
present a comparison of a number of visual and photometric 
properties of galaxies in Figure where the r-band effec¬ 
tive radius r-band ellipticity, r-band absolute magni¬ 
tude Mr, rest-frame and extinction corrected u — r colour, 
and r-band Sersic index rir of galaxies in filaments, tendrils, 
and voids (in blue, green, and red respectively) are com¬ 
pared to each other. The figure is arranged as a grid, with 
the scatterplots showing relationships between these various 
parameters, while the ID histograms at the end of each row 
represent the distribution of that single parameter for the 
three populations of galaxies. 

When we produce our mass normalised samples, we do 
not attempt to normalise for location within groups (i.e., 
dark matter halos), so it is plausible that the fraction of 
satellite galaxies differs between the various samples. How¬ 
ever, as discussed below, this does not appear to produce 
artificial differences between the void, tendril and filament 
samples. 


No population of galaxies uniquely inhabits a particular 
corner of the parameter spaces we are studying (for example, 
filament galaxies are not exclusively red, and do show some 
signal of bimodality in colour). In the context of the proper¬ 
ties shown in this figure, large scale structure has little im¬ 
pact on galaxy properties, beyond determining the form of 
the GSMF which is an important environmental effect in it¬ 
self. There are no statistically significant differences between 
the colour and Sersic index of galaxies in filaments, tendrils, 
and voids; performing a two-sample KS-test between the 
u — r and Ur values between all three populations yields P 
values far exceeding 0.05, indicating that they are consistent 
with having been drawn from the same distribution. We only 
see subtle differences with respect to their effective radius 
and ellipticity. The difference in ellipticity might therefore 
be an indicator of discs being more prevalent in voids. The 
absolute magnitudes of void galaxies suggests that they are 
slightly, but consistently fainter than galaxies in filaments 
and tendrils. Given that this is a mass controlled sample, 
this implies a systematically higher mass-to-light ratio for 
void galaxies and is consistent with the picture of galaxies 
in more dense environments being able to acccumulate a 
larger, brighter stellar population. 


4.3 Total energy output 


The cosmic spectral energy density (CSED, Driver et al. 
2013b I describes the total energy being generated by galax¬ 


ies within a representative volume of the Universe at a given 
epoch. The CSED provides an empirical measurement of the 
energy budget of a population of galaxies by wavelength, 
which is instrumental in understanding the different physics 
and stellar populations of different types of galaxies. An at¬ 
tractive feature of the CSED is that it can be modeled (e.g. 
Somerville et al. 2012 Driver et al. 2013a I given a priori 


knowledge of the cosmic star formation history, initial mass 
function, a stellar population model, and a dust attenuation 


model. Following the prescription of Driver et al. (2013bI, 


we construct CSEDs of galaxies in filaments, tendrils, and 
voids across 20 bands in the FUV-FIR range. We focus on 
this wavelength range as it is entirely dominated either by 
direct starlight (FUV to NIR) or reprocessed starlight that 
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Figure 6. 1 and 2D distributions of galaxies as a function of effective radius, ellipticity, absolute magnitude Mr, u — r colour, and r-band 
Sersic index rir for mass normalised galaxies in filaments, tendrils and voids (shown as blue, green and red lines and points respectively). 
Each frame plots two of these parameters against each other, while the histograms show the distributions of each one individually. Aside 
from very subtle differences in the distribution of ellipticity, effective radius, and absolute magnitude, there is very little difference in the 
properties of galaxies in the 3 environments. 


has been re-emitted by dust (FIR). At 2 : ^ 0.1 the impact 
of AGN at these wavelengths is neglible (Driver et al.|2008 
2013b||. 


We construct our CSEDs by first summing up the lu¬ 
minosities of galaxies with photometry in all 20 bands that 
reside in filaments, tendrils, and voids; we consider only 
galaxies that are observed in all bands. This sum gives the 
luminosity density of those galaxies, which can then be con¬ 
verted to an observed energy density Cobs by multiplying by 
the effective mean frequency of each filter, or band. The in- 
trinsic energy density is then tint = eobs/Peac.A, where Peac,A 
is the mean photon escape fraction for each band, derived 


as in Driver et al. (2007 \ using the dust attenuation model 


of Tuffs et al. (20041 (see also Popescu et al.|2011 1 . Finally, 
we normalise the energy densities in all bands to that of the 
JC-band, and show the results in Figure The data points 
represent the flux in each band, and we use MAGPHYS 
(Multi-wavelength Analysis of Galaxy Physical Properties; 
|da Cunha et al.|2012 1 to fit an unattenuated CSED to each 
data set; these fits are designed to be for illustrative pur¬ 
poses only. 


Figure shows some differences between the energy 
outputs of galaxies in different large scale environments. 
Void galaxies show a very subtle decrease in emission in the 
FIR/submm bands. It is clear from the fits shown in this 
hgure that the dust content is very poorly constrained, and 
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Figure 3. Galaxy stellar mass functions of galaxies in different 
environments by 0.1 dex bins in log stellar mass, with single- 
Schechter function fits to each. The region shaded in grey marks 
our mass selection limit. Each panel shows GSMFs for galaxies in 
different types of environment: large scale structure (top), groups 
(middle), and pairs (bottom), with the total GSMF for all galaxies 
shown as the black line in the top panel. Each set of GSMFs adds 
up to this total set. See Figure 3 in |Kelvin et al.| l |2014[ | for GSMFs 
of galaxies by their morphological type. 


so we cannot make statistically significant statements about 
the meaning of this small difference. 


4.4 Metallicity 


Measuring the abundance of metals in the gas of galaxies is a 
strong indicator of the chemical reprocessing that has taken 
place in their stars, as well as the presence of heavier ele¬ 
ments in their interstellar medium that have accumulated 
as a result of previous generations of star formation. Dif¬ 
ferent elemental abundances can be affected by a number 
of dynamical processes, and are sensitive to environment, 


as well as stellar mass (e.g. Foster et al. 2012 Lara-Lopez 


9.5 10 10.5 11 11.5 



Figure 4. Kernel density estimates of the stellar mass distribu¬ 
tion in mass normalised galaxy populations (colour coded in the 
same way as Figure]^. The bumps in density are due to Poisson 
noise in the void sample, which are then reproduced in the other 
populations by design. 



loQio M. / h'^ 


Figure 5. Fractional distribution of mass-matched galaxies in dif¬ 
ferent environments as a function of stellar mass bins. Note that 
the fractions are calculated within each environment type (for ex¬ 
ample, the fraction of galaxies in pairs is given as A^pair/(-^pair + 
N’notpair))- Errors shown are calculated by sampling the Beta dis¬ 
tribution < |Camer(m |2013| l. 
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Figure 7. Total flux densities in 20 bands for mass normalised 
galaxies in filaments, tendrils, and voids, shown in blue, green and 
red respectively. The CSED shown in grey has been taken from 
[Driver et al.| < |2013b| | and shows an example of a typical GAMA 
galaxy SED. It has been scaled to lie above the points. 


Figure 8. Metallicities and stellar masses from 
EmLinesGAMAIIvOl and StellarMassesvlB shown for mass 
normalised filament, tendril and void galaxies in blue, green and 
red respectively. 


|et al.|[2013[). Recent results show that galaxies in clusters 


have a tendency to be slightly more metal-rich (Cooper et al. 


2008), but that this relationship is largely driven by galax¬ 


ies in very high density regions and local environment as 
opposed to cluster membership (Ellison et al.|2009). On the 


other hand, Hughes et al. (2013) find that in a sample of 


260 nearby late-type galaxies, the stellar mass-metallicity 
relation is insensitive to environment. 

We examine the mass-metallicity relation for star form¬ 
ing galaxies in filaments, tendrils, and voids using gas metal¬ 
licities. While there is no notable difference between filament 
and tendril galaxies, we do find that void galaxies have sub¬ 
tly lower metallicities in comparison to these other two pop¬ 
ulations, as can be see on the right hand panel of Figure]^ 
We suggest therefore that material in filaments and tendrils 
may be being processed in a similar way, and perhaps less 
efficiently in void galaxies. However, considering the similar¬ 
ities in their CSEDs as shown in Figure this effect is not 
thought to be significant. 


4.5 Comparison to group and pair classifications 

It is instructive to repeat some of the above analysis on 
galaxies classified according to their presence within a group, 
or a galaxy-galaxy pair. This allows us to gain some under¬ 
standing into a ‘hierarchy’ of how environment affects galaxy 
evolution: from pair-pair interactions, to the mass of the halo 
or group in which the galaxy resides, to its presence within 
large scale structure. 

We take our original sample of galaxies and classify 
them not by large scale structure, but by presence within a 
group. If the galaxy is in a group, then it is classified accord¬ 


ing to the mass of the group (taken from GroupFindingv07). 
These are dynamical mass estimates of each group calcu¬ 
lated using the group’s velocity dispersion and radius. See 


Robotham et al. (20111 for a discussion on how these proper¬ 


ties were measured. We also identify all galaxies in pairs that 
are within 100 h~^ kpc fro m each other, with velo city sep¬ 
arations of ^ 1000 kms“^ ( Robotham et aL]|2014 |. 50.16% 
of galaxy pairs in groups exist within groups with 3 or more 
members; in other words, they are pairs within larger groups. 
Finally, note that all samples are mass matched to the void 
galaxy population. 

We take the same properties used to generate Fig¬ 
ure and this time plot them for galaxies in groups in 
high mass, mid mass, and low mass halos (defined as 


Mh > 


‘A40, 10^^ ^ Alrj ^ lO^^h ^AIq and Af/r 


10 ^^h~^MQ respectively), galaxies not in groups, central and 
satellite galaxies in groups, galaxies within and outside pairs, 
and the dominant and sub-dominant galaxy of each pair (all 
of these populations are mass controlled, as described in Sec¬ 
tion |4^. Figure [^displays the distribution of these param¬ 
eters as a function of different environments, with each row 
representing a different classification type. From top to bot¬ 
tom, the rows show classifications by large scale structure, 
presence in a group, and presence in a pair. Pair-|- and Pair- 
denote the dominant and sub-dominant member of a pair. 

Figure [^reveals a number of interesting trends: namely 
that different environments tend to separate out certain 
galaxy properties better than others. An obvious example 
of this is that halo mass (or group membership) does not 
appear to be a strong predictor of galaxy luminosity, or ellip- 
ticity. Similarly, halo mass alone does not strongly correlate 
with Sersic index. There is an overall tendency for galaxies 
that are not in pairs, and not in groups, to exhibit similar 
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Figure 9. The probability distribution functions of effective radius, ellipticity, r-band absolute magnitude Mr, u — r colour, and Sersic 
index Ur for galaxies split by three different environment types: large scale structure (top row), halo mass and galaxies not in groups 
(second row), and paired and unpaired galaxies (bottom row). There is an overall trend for galaxies outside of groups, in group centres, 
and outside of pairs to show similar characteristics, particulary in Re, u — r colour, and rir- 


properties: fainter, bluer, higher ellipticity, and a Sersic in¬ 
dex consistent with being a disc galaxy. On the other hand, 
galaxies that reside within groups and pairs have a bimodal 
colour distribution with a higher red peak, higher Sersic in¬ 
dices consistent with early type galaxies, but do not tend to 
be more luminous; though we make this statement whilst ac¬ 
knowledging the caveat that we have controlled our sample 
for mass. 


Our results regarding the properties of pairs are in ex¬ 
cellent agreement with [Robotham et al.| |2013[ ), who con¬ 
clude that the net effect of close pairs is a supression of star 
formation. This is reflected in our results as paired galax¬ 
ies being far redder than unpaired galaxies (which is also 
in good agreement with Davies et ah, in prep). These re¬ 
sults suggest a scenario where the process of infalling into a 
higher density structure generates changes in the morphol- 
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ogy and colour of a galaxy. This is consistent with the re- 


cent results (e.g. 

Valentinuzzi et aL|2011 Lopes et aL|2013 

Robotham et al. 

2013 Lacerna et al. 20141 which indicate 


that local density plays a more important role in galaxy 
evolution than the mass of the halo in which the galaxy re¬ 
sides. Central galaxies in groups occupy the same parameter 
space as ungrouped galaxies; this is consistent with the find¬ 


ings of Robotham et al. (20131. Ungrouped galaxies can be 


interpreted as the central galaxies of undetected groups (i.e. 
GAMA is not sufficiently deep to observe their satellites). 
Notably, Figure 10 of Robotham et al. (20131 establishes 
that central galaxies with M* « 1.25 x Mq are 

more likely to be late type than centrals of a higher mass. 
The distribution of Sersic indices for central galaxies in our 
sample of 9.5 ^ logM^,/h~^M q ^ 11 is consistent with this 
finding, and suggests that central galaxies with even lower 
masses are even more likely to be late type. It is the satel¬ 
lite galaxies in groups that are most directly affected by the 
group environment, as shown in their distribution in Figure 
^ Ag ain, this is consistent with the findings of |Robotham| 
^ al.| ( [2013 1. 

Figure m shows an alternative view of the properties 
shown in previous sections, and for a wider selection of 
environments. Alongside the environment types discussed 
above, for this analysis we also consider visual classifica¬ 
tion of morphology, i.e. E, SOa, Sbc, and Sd galaxies. These 
galaxies were visually classified by three pairs of indepen¬ 
dent observers as being early type or late type, based on the 
presence of a disc or disc-like feature. The sample contains 
approximately 20,000 galaxies, from 0.002 < 2 < 0.1. All en¬ 
vironment and morphological samples are mass matched to 
the void galaxies. The environment types are colour coded 
according to common ‘types’: morphological types are shown 
in blue; pairs in green; groups in red; and large scale struc¬ 
ture in orange. Shaded regions show the spread of the data, 
and the error bars are the standard error about the median. 


4-5.1 Elliptical fraction 

We can similarly look at the fractions of early and late type 
galaxies in these environments, using the visually classified 
morphological galaxy sample. We identify all galaxies in this 
sample that belong to a group, and track the fraction of early 
type galaxies in groups as a function of halo mass. These 
results can be seen in Figure [TTj where the left panel shows, 
in red, the fraction of early type galaxies as a function of 
halo mass. We do not see a correlation between early type 
fraction and group mass, when considering all galaxies (i.e. 
without mass normalisation). 

The right hand panel of Figure 0 displays the total 
fraction of early type galaxies in different environment types. 
These results are consistent with the galaxy types identified 
in Figure in the sense that galaxies in the field have a 
lower early type fraction. Unlike the left panel of the fig¬ 
ure, which shows the elliptical fraction as a function of halo 
mass for all galaxies, the right hand panel shows that mass- 
controlled galaxies in high mass halos have a higher elliptical 
fraction than those in mid and low mass halos. Given that 
the upward trend in the fraction of early type galaxies with 
halo mass is seen only for a narrow subset of mass controlled 
galaxies, we reason that the presence of massive ellipticals 
in halos of all masses drives up the fraction in the left hand 
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Figure 11. Fraction of ellipticals of all galaxies in 
VisualMorphologyv02 that are in groups as a function of 
halo mass shown in red. The fraction of ellipticals does not 
change as a function of group mass. Note that elliptical fractions 
are calculated in bins of halo mass each containing an equal 
number of galaxies. The right hand panel shows the fraction of 
elliptical galaxies in a variety of environments. The error bars 
show 1(T uncertainties in the population fraction, as calculated 
using the QBETA function (|Cameron|2013^ . 


panel. To check if the increase in elliptical fraction in high 
mass groups is simply due to a larger presence of pair galax¬ 
ies in snch systems, we show the fraction of galaxies in pairs 
in our mass normalised sample as a function of halo mass 
in Figure [T^ We note a small, statistically significant down¬ 
ward trend in pair fraction as a function of halo mass. Figure 
[^indicates that the increased fraction of elliptical galaxies 
in more massive halos (for a mass controlled sample) is not 
related to presence in pairs. 


5 HOW DOES MASS NORMALISATION 
AFFECT OUR RESULTS? 


Our analysis and results so far have been restricted to the 
narrow and noisy mass distribution occupied by our sample 
of void galaxies. In this section, we not only examine the 
properties of galaxies without any mass control, but also 
match them to the smoother Baldry et al. (20121 GSMF. 


5.1 Removing mass control 

In order to better understand the effect of controlling stel¬ 
lar mass in our sample, we replicate the results of Figure 
[^for galaxies without mass normalisation. Removing mass 
normalisation has a number of effects, which can be seen in 
Figure An immediately obvious result is that the galaxy 
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Figure 10. Distributions of the same properties shown in Figure]^ but for a greater number of populations. The points plotted are 
medians of that parameter for each population, with shaded regions showing the standard deviation of each parameter, and the error 
bars representing the standard error of the mean (i.e. ST)/sqrt(N — 1). It is for galaxies classified by morphological type that we see the 
greatest variation in properties; all other populations span the same approximate range (as shown by the vertical lines) but do show 
differences in their medians and the widths of their distributions. 


properties shown in this figure show greater trends with en¬ 
vironment (with the exception of those for galaxies classified 
by morphological type), where we continue to see the same 
trends. If the reintroduction of stellar mass variations leads 
to greater trends in the properties of galaxies classified by 
environment and not by morphology, it stands to reason 
that stellar mass is more influential in governing the param¬ 
eter space that galaxies inhabit, cf. large scale environment. 
A secondary conclusion drawn from Figure [T^ is that mor¬ 


phological classifications of galaxies are strong predictors of 
their properties, even when controlling for stellar mass. 

We can also compare the medians of the distributions 
of each property shown in Figure ^[between the mass nor¬ 
malised and non-mass normalised populations. Such a com¬ 
parison quantitatively highlights the effects of mass normal¬ 
isation. In Figure we plot the ratios of the medians of 
the mass normalised and non-mass normalised galaxies for 
each population, for each of the properties discussed in pre- 
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Figure 13. As with Figure but for galaxies whose masses have not been normalised. Most trends are notably more exaggerated in 
this figure, with the exception of the morphologically classified galaxy population, which is a strong indicator for the dominance of stellar 
mass and local density in influencing these parameters. 


vious figures. The void galaxy sample is the same for both 
samples, and so is not shown in this Figure. Different pop¬ 
ulations appear to show a different amount of susceptibility 
to mass normalisation, as evidenced by the larger value of 
the ratio in those properties. Colour and Sersic index appear 
to be the most affected. Central galaxies exhibit the largest 
departure from unity, once again highlighting how the prop¬ 
erties of central group galaxies differ with mass, as shown in 


these results, we study the mass evo¬ 
lution of the parameters examined in Figure [^across envi- 


Robotham et al. (2013 


To further explore 


ronments. We include all galaxies that were used in generat¬ 
ing the GSMFs in Figure without any mass matching in 
any populations. The relationship between stellar mass and 
these parameters (ellipticity, r-band effective radius re, ab¬ 
solute magnitude Mr, u — r colour, and r-band Sersic index 
Ur) is shown in Figure]^ We also include this relationship 
for galaxies classified by their local density (based on the dis¬ 
tance to the 5th nearest neighbour of each galaxy). We split 
this population into 3 bins: containing the lowest Icr of the 
data, between la and 3a, and above 3a. Galaxies in large 
scale structure, groups, and pairs tend to to show a similar 
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Figure 14. Ratios of the median values of each of two galaxy samples shown in Figures |l0| and |l3| (i.e. the ratio between the median 
for the mass normalised sample and the non-mass normalisation sample). The errors shown are standard errors about the median. We 
do not show void galaxies in this Figure, as the population is the same for both samples. 
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Figure 12. Fraction of galaxies that are considered to be in pairs 
as a function of halo mass, in bins of equal size. For each bin we 
plot the fraction Afpair/(Afpair + ^notpair) of all galaxies in halos 
of that mass bin, and include both mass matched and unmatched 
galaxies. We note a slight tendency for the fraction of galaxies 
in pairs to fall in higher mass halos; and that this trend persists 
even when we remove the mass matching. 


evolution in their parameters with stellar mass regardless of 
the environment they reside in. We do see a tendency for 
galaxies in more massive halos, and in pairs, to be system¬ 
atically more red and have a higher Sersic index across the 
whole mass range, but this pattern is not seen for large scale 
structure. Galaxies classified by morphology display consid¬ 
erably different trends. We should not be surprised that this 
is the case for ellipticity, n^, and (although it is notable 
that Ur tends to increase with mass for galaxies classified as 
ellipticals or SOa galaxies). 


9.5 10 10.5 11 11.5 12 12.5 



Figure 16. GSMFs of galaxies separated by environment, as in 
Figure this time compared to the GSMF from [Baldry et al.| 
l|2012|l, which is shown as the dashed magenta line in each panel. 


5.2 Baldry 2012 GSMF 


The Baldry et al. (20121 GSMF describes the mass disti- 


bution of 5210 galaxies from the first three years of the 
GAMA survey, covering an area of 143 square degrees. The 
galaxies are selected such that they are magnitude complete 
(r = 19.4 mag for galaxies in GOO and G15 and r = 19.8 
mag for G12) and with 0.002 < z < 0.06. The GSMF is then 
computed using a density corrected 1 /Fmax method (see Sec¬ 
tion 2.5 of Baldry et al. (20121 for details). We show this 
GSMF as the dashed magenta line in each panel of Figure 
|16[ where the mass functions of galaxies in different environ¬ 
ments, sampled such that they now match to this GSMF are 
also shown. Void galaxies are necessarily removed from this 
analysis. These stellar mass functions have slightly higher 
number densities at higher masses; this is caused by small 
numbers of high mass (^ h~'^M q) galaxies influenc¬ 

ing our Schechter function fits. Across all environments, our 


sample contains only 17 such high mass galaxies, so this 
discrepancy has no discrenible influence on our conclusions. 

The resulting mass matched sample is similar to the 
void matched sample used in previous sections of this paper 
in that it provides a mass-unbiased selection of galaxies. By 


matching to the Baldry et al. (20121 GSMF; however, we are 


able to include a larger number of galaxies at higher mass; 
most importantly, galaxies whose mass is greater than M*. 
It is therefore a high mass alternative to the void matched 
sample used earlier in this paper. A version of Figurej^made 
with this higher mass matched sample is shown in Figure 
|17| Similarly, Figure is equivalent to Figures and 
Large scale structure (i.e. filaments and tendrils) continue 
to have similar characteristics at these higher masses. For 
groups and pairs, the inclusion of higher mass galaxies has 
a number of discernible effects in the distribution of prop¬ 
erties. While group mass continues to be a poor predictor 
of Mr and Ur, centrals, satellites, and ungrouped galaxies 
show a greater variation in their luminosities compared to 
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Figure 15. Relationships between stellar mass and ellipticity, r-band effective radius r^, absolute magnitude Mr, u — r colour, and 
r-band Sersic index rir for galaxies classified according to large scale structure, group membership, pair membership, surface density, and 
morphology. Each population is binned in stellar mass, and the x-coordinate of each bin is the average stellar mass of the galaxies in that 
bin. Errors shown are standard errors about the mean (points with no error bars are single detections). There is an overall tendency for 
galaxies in all environments to be corellated with mass in a similar way, with only morphological classifications showing distinct trends. 


galaxies in groups. All populations show an increase in the 
bimodality of their colours, as a result of the inclusion of 
higher mass galaxies. Central galaxies are now less likely to 
have Sersic indices consistent with having a disk, which is 
consistent with the findings of Robotham et al. (20131, who 
show that higher mass central galaxies are more likely to be 
elliptical. Dominant and sub-dominant pair galaxies show a 
large difference in their Mr distributions compared to Fig¬ 
ure and we see fewer blue dominant pair galaxies. These 
results are also consistent with Robotham et al. (20131, who 
show that more massive dominant pair galaxies are more 
likely to be early type than their lower mass dominant pair 
counterparts. 


6 SUMMARY AND CONCLUSION 

We have used structural and photometric data products 
from the GAMA survey to seek out trends in galaxy proper¬ 
ties as a function of nonlocal environment. We focus in par¬ 
ticular on parameters directly affected by changes in stellar 
populations: colour, brightness, morphology, and total en¬ 
ergy output, as well as some derived parameters. Figure 
highlights the fact that galaxies in different environments 
have strongly varying stellar mass functions. We control for 
this aspect by normalising our samples so that their stel¬ 
lar mass functions match those of void galaxies, which have 
the most extreme mass distribution. This ensures that our 
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Figure 17. As with Figure]^ but for galaxies matched to the|Baldry et al.|||2012[| GSMF. 


comparisons are free from any stellar mass bias. By looking 
at the total CSED of galaxies in different environments in 
Figure we can see that void galaxies have slightly lower 
emission in IR bands, but these fits are poorly constrained. 
Similarly, metallicity does not appear to be strongly affected 
by large scale structure (Figure]^. Figure[^establishes that 


when controlling for the mass distribution, the properties of 
galaxies in different large scale environments are largely the 
same. 

We broaden onr investigation and look at how the prop¬ 
erties examined in Figure [^change for galaxies in group and 
pair environments. Our results are shown in Figure[^and|10[ 
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Figure 18. As with Figures [To| and [T^ but for galaxies matched to the|Baldry et al.|||2012[| GSMF. 


and suggest a tendency for galaxies outside of structures (i.e. 
galaxies not in groups, and not in pairs) to inhabit similar 
parameter spaces. These unpaired and ungrouped galaxies 
are blue, faint, and late type; but we do not see this for 
galaxies in voids. Curiously, group mass does not appear 
to correlate strongly with certain parameters: most notably 
Sersic index, absolute magnitude, and ellipticity. Ungrouped 
galaxies are strongly unimodal in their colour distribution, 
suggesting a change in stellar populations driven by infall 
into groups. We see distributions for central and satellite 
galaxies in groups, as well as dominant and sub-dominant 
galaxies in pairs, that are consistent with other GAMA re¬ 


sults from [Robotham et al.| ( [2013p . The most visible trends 
are for galaxies classified by pair membership, in that galax¬ 
ies in pairs show signs of being much more red (i.e. having 
suppressed star formation), and have a higher Sersic index. 
These trends are seen again in Figure [m which shows the 
fraction of galaxies that are early type to be much higher in 
pairs and higher mass groups (when one controls for stellar 
mass). 


Aware that mass matching to the void galaxy sample 
restricts our results to a narrow mass range, we generate a 
second set of mass controlled galaxies, this time matched to 
the GAMA GSMF from Baldry et al. (20121. While galaxy 
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properties continue to be unaffected by large scale structure 
for this higher mass matched sample, the distributions of 
certain galaxy types reflect the inclusion of higher masses. 
Central galaxies and dominant pair galaxies are now more 
likely to be late type, and all populations show a greater 
bimodality in their colour distributions; these results are 
consistent with the findings of Robotham et al. (20131. We 
further explore the effects of controlling for mass by remov¬ 
ing it entirely for Figure where it is shown that galaxies 
in environments have a tendency to be correllated in similar 
ways as they go to higher stellar masses. The most mas¬ 
sive galaxies in the Universe look the same regardless of 
where they are. Furthermore, removing mass normalisation 
also enhances trends in galaxy properties with structure as 
in Figure [T^ indicating strongly that stellar mass (or what 
drives stellar mass) is the dominant force in galaxy evolu¬ 
tion. We see the most significant effects when looking at pair 
classifications, and is consistent with findings in |Robotham| 


et al. (20131 and Robotham et al. (20141. A macroscopic 


view of our results suggests that the properties of galaxies 
reviewed in this work correllate most strongly with stellar 
mass. 

The results presented in this work serve as a broad 
overview of how large scale structure and environment di¬ 
rectly impact on galaxy evolution, and serve as a broad 
overview of the trends of galaxies in different environments. 
These results motivate further study, in particular to de¬ 
tect mass flows within filaments by analysing the properties 
of spiral galaxies. This field stands to gain a much deeper 
understanding of the environmental mechanisms that influ¬ 
ence galaxy evolution as more detailed surveys of gas in 
more distant galaxies become available with the advent of 
surveys such as ASKAP DINGO, followed by the SKA. The 
GAMA Large Scale Structure Catalogue presents a unique 
opportunity to study these processes within a well defined 
sample of galaxies in extremely underdense environments. 
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